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• Intersubband transitions in a superlattice under homogeneous electric field is studied within the 

tight-binding approximation. Since the levels are equi-populated, the non-zero response appears 
beyond the Born approximation. Calculations are performed in the resonant approximation with 
scattering processes exactly taken into account. The absorption coefficient is equal zero for the reso- 
nant excitation while a negative absorption (gain without inversion) takes place below the resonance. 
A detectable gain in the THz spectral region is obtained for the low-doped GaAs-based superlattice 
(— I ' and spectral dependencies are analyzed taking into account the interplay between homogeneous and 

Ch . inhomogeneous mechanisms of broadening. 
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The examination of stimulated emission due to intersubband transitions of electrons (monopolar laser effect), which 
^ has been carried out during the previous decade, have resulted in mid-IR lasers (see Refs. in Recently, the 

THz laser has also been demonstrated [3|, Ij, |5|, |^ . The standard laser scheme based on vertical transport through 
Q ' the quantum cascade structures, which incorporates the injector and active regions, has been used in both cases. 
Population inversion appears in the active regions and leads to stimulated emission for the mode propagating along 
mid-IR or THz waveguide. In contrast to this, the vertical current in a biased superlattice (BSL) with the Wannie- 
Stark ladder, which appears under the condition 2T ^ es,B 7] (here Ebb/^ is the Bloch frequency and T stands for 
the tunneling matrix element between adjacent QWs), does not change the populations of the levels. Due to this. 



the consideration based on the golden rule approach gives a zero absorption. At the same time, for the wide minigap 

O T~r 

SL, with the width 2T ^ egs, a negative differential conductivity, i.e. gain due to Bloch oscillations, takes place 

o 



This contradiction and the question about THz gain without inversion are discussed in Ref. Last year, agreement 
' between the numerical results for the wide-miniband and hopping regimes of high-frequency response was noted in 

S: 0. 

'"^ ' Since there is no well-defined dispersion relation between energy and momentum, e and p, beyond the Born 
Q ' approximation, one has to consider the intersubband transitions based on the spectral density function, A^{-p), which 
is a finite- width peak [ll|. Let us consider first the two- level model with an identical distribution finction for both 
^ • levels, /e. We take into account the off-resonant transitions with a non-zero detuning energy Ae = fiuj — e^b with 
■ respect to the level splitting energy, s^b, see Fig.l. The intersubband absorption is given by 



(1) 



moreover the relation oc sign(Ae) is obtained under the replacement e — Ae e. Taking into account that 
decreases with e, one immediately obtains uas < if Ae < |l2| . i.e. a gain appears in the BSL with disorder 
beyond the Born approximation. In the Born approximation, when the spectral function is replaced by 5-function, 
one obtains a^e = 0. 

In this paper, we evaluate Eq.(l) for the low-doped BSL taking into account the scattering processes exactly. 
Numerical estimates are performed below taking into account the interplay between homogeneous and inhomogeneous 
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FIG. 1: Off-resonant intersubband transitions (a) and corresponding spectral density functions (6). The dashed curves show 
distribution functions and arrows indicate the transitions with positive and negative detuning energies. 



hrr' = ( ^ 1- Vrx + TEnB ] 5^' + T{5rr'-1 + <^rr' + l), (2) 



mechanisms of broadening. 

Within the framework of the tight-binding approach we describe the electron states in BSL using the matrix 
Hamiltonian: 

2m 

where /2ra is the in-plane kinetic energy operator, m is the effective mass, T/rx is a random potential energy in the 
r-th QW, r = 0, ±1, . . ., which are statistically independent in each QW. The Bloch energy, ebb^ appears in 

(2) due to the shift of levels in the SL with period Z under a homogeneous electric field F, see [l3|- The perturbation 
operator due to a high-frequency transverse field [-Ebj. exp(— iwt) -f c.c] is written as [5hrr' exp{—iu!t) + H.c], where 
the non-diagonal matrix 6hrr> is given by 

dhrr' = —Eiij_{Srr'-l +^rr' + l) (3) 

LU 

and Vfi± = TZ/h. The high-frequency current induced by the perturbation (3), [/^^ exp(— icji) + c.c], is determined 
by the standard formula: 

where 2 is due to spin, spg|| ... is the averaging over in-plane motion, ((. . .)) is the averaging over random potentials 
Vrx, and is the normalization volume. 

The high-frequency contribution to the density matrix in Eq.(4), [6p„, exp{—iujt) + H.c], is governed by the 
linearized equation: 

T ^ 

+ ^{SPr+lr' + SPr-lr' ~ ^Prr'-l " ^Prr' + l) 

+ — Shrr'(pr' - Pr) — 0. (5) 

n 

Here hj. — /2m + Vrx + tsub describes an in-plane motion in the r-th QW and we use the steady state density 
matrix {po)rr' — Srr'Pr i 1-6. wc havc neglected a weak non-diagonal term which is responsible for the tunneling 
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current through the BSL. We restrict ourselves to the consideration of cx contributions only, so that we can omit 
oc T addendums in Eq.(5). Thus, an independent equation for Sp,. = Sp^^i^ takes the form: 

— lOJOp^ + -(hr±lOPr ~op^ hr) 

— -i^-^^Ee,±{pr±i - pr). (6) 

Note that for the coUisionless case {pr±i — Pr) — > 0, so that the response vanishes and is only non-zero due to 

differences in scattering processes for adjacent QWs. In addition, for the resonant approximation, Ihui — ehbI ^ £bs, 

-fl(-) 

one can neglect the contribution of Sp^. 

Writing spgn ... in Eq.(4) in the coordinate representation, we obtain/^ = {i2evs±/L^){J2rI dx(5p'^(+^(x,x)). Next, 
we describe the electron states in the r-th QW by the use of the eigenstate problem {p^ /2m + K-x)V'rx ~ ^ri/V'rxj 
where a quantum number v marks an in-plane state. Using this basis, we transform /^^ into 

L3 

and the linearized kinetic equation takes the form 



(7) 



= - /e„„,] j dx^+lxV'.^x *• (8) 

Here A — ^ +0 and we use the quasi-equilibrium distribution pr — fp2/2m+Vr^, where is the Fermi function with 
identical chemical potentials, p, and temperatures, Tg, for any QW. We introduce the conductivity, a^j, according to 
Ilu = o'ujEp,j_, and Eqs. (7,8) give us: 

^ _ ^ 2{ev,^)^ //y iU^,^ - fe,^^,)Q7+l.r \ \ 



where Q'^'^^, = J dx.ip'/.^ip'/.i^ is the overlap factor. Thus, we have evaluated the expression for the response with the 
scattering processes exactly taken into account. 

Below we consider the absorption coefficient introduced according to a^j = (47r/c-\/e)i?e(7tj, where e is the dielectric 
permittivity, which is supposed uniform across the structure. In order to perform avera ging in Eq.(9), we use the 
spectral density function in the r-th QW determined as Are (x, x') = 'ip'^^ip'/.^5{erv — e) JJj , so that aAe is written 
as follows: 

OlAe = T-pr- / de{fe-Ae - fe) (10) 



>^ J dy^J dx'^((A+le(x,x')Ae-Ae(x',x))) 

with huj ~ £6S in the resonant approximation. 

We turn now to averaging over short-range and large-scale potentials taking into account that we are considering 
SL under a homogeneous bias voltage. Due to this the averaged characteristics of scattering processes, both for 
homogeneous and inhomogeneous mechanisms, do not dependent on the QW number r. It is convenient to use the 
Wigner representation and the average of the spectral functions in (10) takes the form: 

dxrfx' 



L3 



^((A+1.(X, x')Ae-As(x', x))) 
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FIG. 2: Dimensionless gain —a/^^/ao versus detuning energy — Ae/7 for the homogeneous broadening case (F = 0). Panels 
(a) and (6) correspond to degenerate (/i/7 =3) and non-degenerate {jJ./'y = -1) electrons. Solid, dashed, and dotted curves 
correspond to Te/'y =0.3, 1, and 3 respectively. 
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TABLE I: 3D concentrations, measured in cm versus dimensionless /i and Te for the cases plotted in Figs. 2 (F = 0), 3 
(7 = 0), and 4 (F = 7). 



Here we took into account that ((. . .)) does not depend on x and L~^^^, = . Performing the averaging over 
short-range potential, we write the spectral function {Are{v^^)) ~ IinG^^f'{p,x.)/TT through the retarded Green's 
function given by: 

G^f(p,x) = (ep-«;,x-e-5])-i. (12) 

Here Wrx is a large-scale part of potential in the r-th QW and S is the self-energy function arising from the short-range 
scattering (see similar calculations in ^3). Below we consider the case of scattering by zero-radius centers when /mS 
does not depend on e, p or x. _ReE, which is logarithmically divergent without a small-distance cutoff, is included into 
the detuning energy Ae, so that the only homogeneous broadening contribution, —i"f, appears in the denominator of 
Rq. (12). Performing the averaging over large-scale potentials we write the spectral density in the integral form: 

A,(ep) ^ y° ^eKe,-e-^',)t/n^-i^t/nf/2 ^ ^^^^^ (^3) 

where F = \/ (w^rx) is the inhomogeneous broadening energy. Using the in-plane isotropy of the problem, we finally 
transform Eq.(lO) into 

«Ae^^ r- / de (14) 

he ^ee<i,BZ J Jo 

X^e-Ae/2(0^e+A£/2(C)(/e-Ae/2 " /e+Ae/2)- 

Note, that for the coUisionless case the product of spectral functions under the integrals is transformed into — e + 
Ae/2)5(e-e- Ae/2). 

Further, we calculate the spectral dependencies aa^elc^o given by Eq.(14) with the use of the quasi-equilibrium 
Fermi distribution /e. The characteristic absorption ao is introduced here according to ao — [e^ l'hc)mw\^l ^e\mZ 
and = — ck_Ae, so that we consider only the region Ae < 0. First, we examine the cases of homogenous (F = 0) 
and inhomogeneous (7 = 0) broadening. The dimensionless gain is plotted for the cases of degenerate [Figs. 2(a) 
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FIG. 4: Dimensionless gain —a/^n/o-o versus — Ae/2(7 + r) for the case 7 = F. Panels (a) and (6) correspond to degenerate 
{lJ./2"f =3) and non-degenerate {n/2'y = -1) electrons. Solid, dashed, and dotted curves correspond to re/27 =0.3, 1, and 3 
respectively. 



and 3(a)] and nondegenerate [Figs.2(&) and 3(6)] electrons with the concentrations given in Table I. According to 
these data, concentration increases with Te for all cases but the peak gain varies in a different manner: quenching 
or enhancement of gain occures for degenerate or non-degenerate electrons respectively. High-energy tails of gain 
and lower peak values take place for the homogeneous broadening case, when (13) has a Lorentzian shape. For the 
inhomogcncous broadening case, Ae(^) has a Gaussian shape and appears to be a sharper function. In Fig. 4 we 
present the case 7 = F, plotting the dimensionless gain versus A£/2(7-|-r), where 2(7 + F) is the total width of (13). 
One can see both the same style of spectral dependencies and the same temperature/concentration dependencies. 

Next we turn to estimates of the maximal gain for GaAs/AZo.sGao.y^s-based BSL with T =0.5 meV, corresponding 
to the barrier width of 6 nm and Z =15 nm. For the level splitting energy e%B = 10 meV, which is correspondent 
to the transverse field F =6.7 kV/cm, one obtains ao ~6.6 cm~^, so that the peak gain appears to be between 5 
and 20 cm~^ for different parameters used in Figs. 2-4. Note, that ao oc in the framework of the tight-binding 
approximation and gain increases rapidly for the narrow barrier case, for example gain exceeds the experimental data 
if T =1 meV. 

Let us discuss the main assumptions used. The tight-binding approach is valid under the condition eg^ 3> IT 
which is satisfied for the numerical estimates performed; note, that beyond the Born approximation the broadening 
can be comparable with the electron energy determined through /i and Tg. We restrict ourselves to the case of homo- 
geneous field and concentration distributions neglecting a possible domain formation due to the negative differential 
conductivity at low frequencies [isl ]. One can avoid instabilities in a short enough BSL because the THz modes 
propagate in the in-plane directions. In spite of the general expressions (10-12) are written through an arbitrary 
self-energy function E, the final calculations were performed for the model included scattering by zero-radius centers 
and large-scale potential. Such a model describes the interplay between homogeneous and inhomogeneous broadening 
with the use of statistically independent random potentials in each QW. The Coulomb correlations, which modify the 
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response as the concentration increases, are not taken into account here. This contribution, as well as consideration 
of intermediate-scale potential, require a special consideration in analogy with the case of a single QW . 

In conclusion, we have considered the resonant intersubband response of a BSL and have described gain without 
inversion beyond the Born approximation. It seems likely that this contribution can be found experimentally and 
more detailed numerical calculations are necessary in order to estimate a potential for applications. 
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